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We report on a systematic study of the superconducting characteristics and Pr crystalline- 
electric-field (CEF) levels of the filled skutterudite PrRu 4 Asi 2 (T c = 2.33 K). The temperature 
dependences of the upper critical field H C 2 and the Ginzburg-Landau (Maki) parameter K2 
suggest s-wave clean-limit superconductivity. The electronic specific heat coefficient 7 ~ 95 
mJ/K 2 mol, which is ~ 1.5 times larger than that of LaPauAsu, indicates 4/-originating quasi- 
particle mass enhancement. The magnetic susceptibility x(T) indicates that the CEF ground 
state is a T\ singlet and a F^ 1 ' triplet first excited state lies at Acef ~30 K above. A system- 
atic comparison among PrOs4Sbi2, PrRu4Sbi2, PrRu4Asi2, and La-based reference compounds 
suggests that the inelastic exchange scattering and aspherical charge scattering of conduction 
electrons from CEF-split 4/ levels play an essential role in the quasiparticle mass enhancement 
and in determining the value of T c in the Pr-based filled skutterudites. 

KEYWORDS: Filled skutterudite, PrRu4Asi2, specific heat, susceptibility, superconductivity, crystal electric 
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Among heavy-fermion superconductors, the filled skut- 
terudite PrOs4St>i2 has attracted considerable atten- 
tion for its unconventional superconducting (SC) prop- 
erties. 1,2 Some of the characteristic features are bro- 
ken time-reversal symmetry, 3 odd parity, 4 possible struc- 
ture in the SC gap, 5,6 multi band superconductivity, 7 
and adjacent quadrupole ordering and its fluctuations 
(quadrupole excitons). 8-11 In contrast, PrRu 4 Sbi2 is a 
conventional BCS-type superconductor as suggested by 
a coherence peak in the NQR spin-lattice relaxation rate 
l/Ti. 12 Systematic understanding for such variety in the 
SC properties of Pr-based filled skutterudites is one of 
the important issues (see Table II for a summary of the 
SC parameters). 

It has been pointed out that the crystalline electric 
field (CEF) level scheme of Pr ions may play an essential 
role. In the Th site symmetry, the J = 4 multiplet of Pr 3+ 
ions splits into four sublevels, 13 namely, a singlet r 1; a 
non-Kramers nonmagnetic doublet T23, and two triplets 

and 1^ . PrOs4Sbi2 has Ti ground state with a 
first excited state with the energy separation Acef = 8 
K. 11,14 In contrast, PrRu 4 Sbi 2 has Ti - levels with 
Acef = 65 K. 15-17 The differences in the type of the 
first excited state and Acef may be a key reason for the 
different SC properties. 

In this paper, we report the specific heat and magnetic 
susceptibility measurements of another Pr-based filled 
skutterudite PrRu4Asi2 18 to study the SC properties and 
the CEF level scheme. The results indicate that the su- 
perconductivity is of the conventional BCS type with a 
moderate quasiparticle mass enhancement. A systematic 
comparison among these Pr-based compounds and their 
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Table I. Fractional coordinates and Debye-Waller factors U cq of 
PrRu4Asi2 determined from the single-crystal x-ray diffraction 
data. 



atom 


site 
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C/oq(A') 


Pr 


2a 











1.082(4) 


Ru 


8c 


1/4 


1/4 


1/4 


0.294(3) 


As 


24g 





0.14954(3) 


0.35019(3) 


0.387(5) 



La-based references reveals that inelastic scatterings of 
conduction electrons from CEF-split 4/ levels play an es- 
sential role for the SC properties and quasiparticle mass 
enhancement in the Pr-based filled skutterudites. 

Polycrystalline PrRu4Asi2 samples are synthesized us- 
ing stoichiometric amounts of 3N(99.9% pure)-Pr, 4N- 
Ru, and 5N-As powders by the high-pressure and high- 
temperature method (~ 4 GPa at 900 °C). 19 Powder x- 
ray diffraction indicates the inclusion of a small amount 
of the RuAs2 phase; the maximum peak height of RuAs2 
is approximately 10 % of that of PrRu4Asi2- The struc- 
tural parameters of PrRu4Asi2 are determined by single- 
crystal x-ray diffraction using an imaging plate detec- 
tor with Mo K Q radiation. A small single crystal with 
approximate dimensions of 0.10 x 0.07 x 0.05 mm is 
selected from the ingot. The structural parameters are 
successfully refined based on the 179 independent reflec- 



/E\Fo\ 



tions with the agreement factor ^2 

= 0.023, where F and F c are the experimental and cal- 
culated structure factors, respectively. Fractional coordi- 
nates and Debye-Waller factors {U eq ) are listed in Ta- 
ble I. The obtained lattice constant of 8.491(4) A agrees 
with a reported value. 20 t/ eq (Pr) is ~ 3 times larger than 
Z7 eq (Ru) and U cq (As). Note that Z7 eq (Pr) in PrRu 4 As 12 
is larger than that in PrRu4Pi2, 21 indicating larger ther- 
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Fig. 1. (Color online) Specific heat divided by temperature (C/T) 
of PrRu4Asi2 measured in several magnetic fields. The bro- 
ken curve represents the phonon part C p h/T determined for 
LaRu4Asi2. 18 



mal vibrations of Pr ions in Asi2 cages, whose size is 
larger than that of P12 cages. 

The specific heat C (H, T) is measured by a quasiadi- 
abatic heat pulse method 8 using a dilution refrigerator 
equipped with an 8 T SC magnet. The temperature in- 
crement caused by each heat pulse is controlled to ~2 
% for the usual measurement and to ~0.5 % in limited 
temperature ranges where the SC transition occurs. The 
magnetic susceptibility x is measured in 1.8 < T < 300 
K using an SC quantum interference device (SQUID) 
magnetometer (Quantum Design Inc.). 

The temperature dependence of C/T measured in se- 
lected magnetic fields is shown in Fig. 1. In zero field, a 
clear jump appears at the SC transition temperature T c 
= 2.33 K, which is very close to a resistively determined 
value of 2.4 K. 18 With increasing field, T c shifts to lower 
temperatures and the size of the jump becomes smaller 
(for details, see the inset of Fig. 2). 
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Fig. 2. (Color online) Upper critical field /UnHc2 vs T determined 
by T-sweep and H-sweep measurements of C. The error bars 
denote the 10% - 90% width of the specific heat jump. The 
broken line shows the best-fit WHH curve (see text). The inset 
shows C ei /T vs T near T c . 



In the measured temperature range, C consists of three 
terms, namely, the electronic contribution (C e i), phonon 
contribution (Cph), and nuclear Schottky contribution 
(C n ). In zero field, C n = A(H)/T 2 appears as a low- 
T upturn below 0.7 K. This term gradually develops 
with increasing field. The phonon term C p h = /3T 3 with 



= 1.53 mJ/K 4 mol determined for LaRu4Asi2 18 is indi- 
cated by the broken curve in Fig. 1. By the fitting using 
C(T, H) = C e i(T, H) + A(H)/T 2 + /3T 3 , C el (T, H) is ex- 
tracted (see the inset of Fig. 2). 

Note that there appears a slight jump at Ta ~ 6 K in 
zero field. This anomaly is probably due to a magnetic or- 
dering in a secondary phase included in the present sam- 
ple since no corresponding anomalies appear in recent 
75 As-NQR measurements. 22 This phase has not been 
identified by the x-ray diffraction measurements. To syn- 
thesize a sample free from the secondary phase, sample 
preparation is conducted several times under different 
conditions and the best sample is used for the present 
measurements. Probably due to the contamination of the 
secondary phase, C/T decreases with increasing field in 
the normal state, making it difficult to estimate the elec- 
tronic specific heat coefficient (7). Nevertheless, from the 
C data at 8 T, where the contribution of the secondary 
phase is largely suppressed, 7 ~ 95 mJ/K 2 mol is esti- 
mated. The specific heat jump AC/jT c ~ 0.83, smaller 
than 1.43 for a weak-coupling BCS superconductor, and 
apparent ~ T 2 dependence in C(T) below T c may sug- 
gest SC gap anisotropy or multi band superconductivity. 
However, a clear coherence peak observed in As-NQR 
1/Ti indicates that such nonuniformity would not be very 
large. 22 

From the C c \ (T, H) data, the upper critical field 
l^oH c2 {T) is determined, as shown in the H-vs-T phase 
diagram of Fig. 2. It appears that the n$H C 2{T) curve is 
well reproduced by the Werthamer-Hefland-Hohemberg 
(WHH) formula for the clean limit 23,24 with the initial 
slope [-d(ji H c2 )/dT] Tc = 0.365 T/K. This fact indi- 
cates that fioH C 2 is mainly determined by the orbital ef- 
fect and the spin Pauli paramagnetic effect does not play 
a dominant role; the Pauli limiting field is /irj-ffp = 3.9 T 
<C /J.oH C 2- The T — value of /io-ff C 2 is estimated to be 
- 0.62 T by the WHH fitting. 

There appears a slight upward deviation in the 
HoH c2 (T) curve below 0.05 T, leading to an enhance- 
ment in T c by 2% in zero field compared with the best- 
fit WHH curve for ^qH > 0.05 T. A similar behavior is 
observed in the case of PrOs 4 Sbi 2 28 and YNi 2 B 2 C, 29 in- 
dicating that the upward deviation in T c may come from 
the multi band effect. 

By combining the fioH C 2 and AC /T c data, the 
Ginzburg-Landau parameter (or Maki parameter) 30 k 2 
is determined using the thermodynamical relation: 



AC/T C 



d((ioH c2 ) 
dT c , 



1 



4tt (24 - 1) p A ' 



(1) 



where (3a is 1.16 for a triangular vortex lattice. 31 The 
result is shown in Fig. 3. The value of k 2 is ~ 10 at 
T = T c (in zero field) and increases with decreasing tem- 
perature. 32 The behavior of dn^/dT < is again con- 
sistent with the negligible Pauli paramagnetic depair- 
ing effect mentioned above; for superconductors with a 
strong Pauli paramagnetic effect, dn 2 /dT > appears 
(e.g., CeCoIn 5 ; see Ref. 33). The temperature depen- 
dence of K2 is affected by impurity scatterings. 34 In the 
clean limit, the theoretical parameter k 2 (T) diverges as 
T -> with k 2 oc y/\ii{T c /T) 30 > 35 as shown in Fig. 3. 
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Table II. Superconducting parameters of Pr-based filled skutterudites and their isostructural La-based reference compounds. 
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First excited state 
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m*/mo for 7-branch 


Reference 




(K) 


(T) 


(T/K) 


(K) 




(m.T/K 2 mol) 






PrRu4Asi2 


2.33 


0.62 


0.37 


~ 30 


r (l) 

1 4 


95 




This work 


LaRu4Asi2 


10.3 


0.72 


0.08 






73 




18 


PrRu4Sbi2 


1.03 


0.2 


0.24 


65 


r (l) 
1 4 


59 


1.6 


15, 16,25 


LaRu4Sbi2 


3.58 


0.28 


0.12 






37 


1.4 


15, 16 


PrOs 4 Sbi2 


1.81 


2.3 


1.9 


8 


r (S) 
1 4 


~ 500 


7.6 


1,11,26 


LaOs4Sbi2 


0.74 


0.04 


0.095 






55 


2.8 


26,27 



clean limit g/L = 0) 



PrRu 4 As 12 




1 1.5 

T(K) 



Fig. 3. Temperature dependence of Ginzburg-Landau parameter 
(Maki parameter) K2 for PrRu4Asi2 . The broken curve indicates 
the theoretical curve for the clean limit. The solid curve is a guide 
to the eye. 



However, in the dirty limit, the temperature dependence 
of K2 is largely suppressed (k 2 (0)/k2(?c) ~ 1.2). The 
K2 (T) data in Fig. 3 indicates that the present sample is 
rather close to the clean limit. 

The magnetic susceptibility x(T) is shown in Fig. 4. 
At high temperatures, x(^) can be well described by 




Fig. 4. Temperature dependences of susceptibility \ (open circle) 
and inverse susceptibility x _1 (open square). The dashed line 
shows the best- fit CEF model calculation (see text). 



a Curie-Weiss law \ — NA^ 2 e ff/ikB{T — Qcw), where 
the effective magnetic moment [i e ff — 3.30/iB/f.u. and 
the Curie- Weiss temperature &cw — -11 K. The ob- 
tained fx e ff value is close to the Pr 3+ free-ion value (3.58 
/is/Pr), indicating a well-localized character of Pr 4/ 
electrons. Below ~10 K, x(T) shows a deviation from the 
Curie- Weiss law, exhibiting a saturation tendency. This 
behavior suggests that the low-lying CEF level scheme of 
Pr 3+ ions has a Ti singlet ground state and a triplet 



(2) 

first excited state. Note that in the case of a r\ ' first ex- 
cited state, as realized in PrOs4Sbi2, a maximum should 
appear in xiT)- 1 ' 8 This interpretation agrees with the 
electronic entropy S C \(T, H) (not shown) calculated from 
the C c i(T, H) data. S el at 8 K is less than 20 % of R In 2, 
ruling out the possibility of any degenerate ground states. 
In the measured temperature range, dS c \/dH is always 
negative. In PrOs 4 Sbi2, dS e \/dH > is observed below 
3 K. 8 This behavior is caused by the Zeeman splitting of 
the low-lying triplet first excited state at Acef = 8 K. 
This fact indicates that Acef ~S> 8 K for PrRu4Asi2- For 
a rough estimation of the low-lying CEF level scheme, a 
least-squares fit has been made based on the Oh Lea- 
Leask-Wolf scheme 36 (omitting the 0\ — 0\ term char- 
acteristic of the Th site symmetry in the CEF Hamilto- 
nian). The best overall fit of the x{T) data (shown by the 
dashed line in Fig. 4) corresponds to a T% ground state 
with a F 4 (r^ in Th) triplet excited state at Acef ~ 30 
K. 37 For an accurate determination of the overall CEF 
level scheme, inelastic neutron scattering measurements 
would be necessary. 

Only for LaOs 4 Sbi2, among the three systems shown 
in table II, an enhancement in T c is caused by the 
replacement of La by Pr ions. This behavior may be 
accounted for by considering the predominant type of 
conduction-electron scatterings from Pr 4/ electrons. 
Fulde et al. 38 have theoretically demonstrated that in- 
elastic aspherical charge scattering (ACS) associated 
with 4/ quadrupole moments leads to an increase in T c , 
while exchange scattering (ES) decreases T c . For a qual- 
itative discussion, a comparison with the theoretical cal- 
culation is made in Fig. 5. For LaOs4Sbi2, since the first 

(2) 

excited state of the replaced Pr ions is T 4 (~ T5 in 
Oh having large off-diagonal quadrupole moments with 
Ti), the ACS may dominate over the ES. Furthermore, 
AcEF/2T c (LaOs4Sbi2) ~ 5 happens to be the best con- 
dition for the enhancement in T c . On the other hand, 
in LaRu4Asi2 and LaRu4Sbi2, the first excited state of 
the replaced Pr ions is T^ 1 ' (~ T4 in Oh having large 
off-diagonal dipole moments with Ti) and the resulting 
predominant ES will cause a strong pair-breaking effect, 
thereby decreasing T c . 

The quasiparticle mass enhancement is also expected 
to be caused by the interactions of conduction electrons 
with low-lying CEF levels. Flude et al. have proposed 
a model for the mass enhancement due to virtual CEF 
excitations: 39 



1 



to 



(gj - iyJi f N(0) 



2|<i|J|j>f 



*CEF 
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Fig. 5. (Color online) Two curves showing calculated changes in 
T c vs Acef/2T c by doping impurity ions with singlet-singlet 
CEF levels (see Ref. 38 for the parameter definition). For a rough 
estimation of the Pr doping effect on T c , Acef/2T c for LaT4Xi2 
is indicated. The expected dominant scattering process is indi- 
cated by solid circles. 



where gj is the Lande factor, J s f is the exchange inte- 
gral coupling the conduction electrons to the / electrons, 
N(0) is the bare conduction electron density of states at 
the Fermi level, and (i \J\j) is the magnetic dipole ma- 
trix element calculated using the derived CEF parame- 
ters. Figure 6 shows several physical quantities that pro- 
vide measures of the mass enhancement as a function 
of Aqef- There appears a rough trend of m*/m — 1 oc 
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Fig. 6. (Color online) Mass enhancement factor m* /mo — 1, es- 
timated from several physical quantities, is plotted as a function 
of A CEF . 



1/Acef in Fig. 6. However, as discussed above for the 
effect on T c , not only the ES process but also the ACS 
process should contribute to the mass enhancement with 
a different prefactor of 1/Acef- In PrOs4Sbi2, the im- 
portance of quadrupolar degrees of freedom is inferred 
from a strong correlation between quadrupolar excitons 
and superconductivity. 11 Further experimental and the- 
oretical studies on the role of quadrupoles in Pr-based 
filled skutterudites will be needed. 
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